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Abstract. Using the technique of the elliptic Frobenius determinant, we construct new 
elliptic solutions of the QD-algorithm. These solutions can be interpreted as elliptic solutions 
of the discrete-time Toda chain as well. As a by-product, we obtain new explicit orthogonal 
and biorthogonal polynomials in terms of the elliptic hypergeometric function :iE2{z). Their 
1^ ' recurrence coefficients are expressed in terms of the elliptic functions. In the degenerate 

\ case we obtain the Krall-Jacobi polynomials and their biorthogonal analogs. 
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0\ ■ 1 Introduction 

O 

, In this paper we present new explicit solutions for the two-point QD-algorithm [5] (which is 

equivalent to the discrete-time relativistic Toda chain |23t l28l I16j). These solutions can be 
naturally constructed starting from the famous Frobenius elliptic determinant (see, e.g., [HI [3]). 
This approach allows one to find an explicit expression for corresponding Laurent biorthogonal 
' polynomials in terms of the elliptic hypergeometric function ■iE2{z). These polynomials contain 

several free parameters and appear to be biorthogonal on the unit circle with respect to a dense 
point spectrum. In two special cases we already obtained explicit examples of cn- and dn- 
elliptic polynomials which are orthogonal on the unit circle with respect to a positive dense 
point measure [33]. These polynomials provide first known explicit (i.e. expressed in terms 
of the elliptic hypergeometric function) examples of such measures (see also [24] for general 
properties of polynomials orthogonal with respect to measures of singular type and [21] for an 
example of such polynomials). The obtained polynomials Pn{z) possess a remarkable "classical" 
property. This means that VPn{z) = UnPn-iiz), where P is a generalized derivative operator: 
T>z^ = iinz"'~^ (with some coefficients Hn) and Pn{z) are polynomials of the same type but 
with shifted parameters. In our case the operator T> is an elliptic generalization of the ordinary 
derivative operator V = dz, Hn = ^ and (^-derivative operator with = (g" — l)/{q — 1). 



*This paper is a contribution to the Proceedings of the Workshop "Elliptic Integrable Systems, Isomonodromy 
Problems, and Hypergeometric Functions" (July 21-25, 2008, MPIM, Bonn, Germany). The full collection is 
available at |http: / /www. emis.de/journals / SIGMA /Elliptic-Integrable-Systems.html , 
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In the degenerated case (when both periods of elhptic functions become infinity) we obtain 
biorthogonal analogs of the Krall-Jacobi orthogonal polynomials. We show that these biorthogo- 
nal polynomials satisfy a 4th order differential equation which can be presented in the form of 
quadratic operator pencil. 



2 Laurent biorthogonal polynomials and their basic properties 

The Laurent biorthogonal polynomials LBP Pn{z) appeared in problems connected with the 
two-points Pade approximations (see, e.g., [H]). 

We shall recall their definition and general properties (see, e.g., [HI [121 , where equivalent 
Laurent orthogonal functions are considered) . 

Let C be some linear functional defined on all possible monomials by the moments 



n = 0,±l,±2.... 



In general the moments c„ are arbitrary complex numbers. The functional C is defined on 

N2 

the space of Laurent polynomials Viz) = o-nZ^ where On are arbitrary complex numbers 

n=-Ni 

and arbitrary integers: 
C{V{z)} = ^ a„c„. 

n=-Ni 

The monic LBP Pn{z) are defined by the determinant [12] 



Pn{z) 
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Co 
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(2.1) 



where A„ is defined as the Toeplitz determinant 



A. 



Co 

C-l 



Cl 
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Cn-1 
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Cl-n C2-n • • • Cq 

It is obvious from definition (j2.ip that the polynomials Pn{z) satisfy the orthogonality property 

C{Pn{z)z-^} = hn6kn, < k < U, 

where the normalization constants hn are 

/iQ = Co, hn = An+l/An- 

This orthogonality property can be rewritten as the biorthogonal relation |22l I12j . 

/:{Pn{z)Qmil/z)} = hn5nm, 

where the polynomials Quiz) are defined by the formula 



Qn{z) = (A„)-i 



Co 

Cl 



C-l 
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Cn-1 Cn-2 
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Cl-n 
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(2.2) 
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We note that the polynomials Qn{z) are again LBP with moments c\t^^ = c_ 
In what follows we will assume that 

A„y^O, n = l,2,... 

and that 

AW^O, n=l,2,..., 

where by An we denote the determinants 



(2.3) 



(2.4) 



A 



1 A^-'') 



Cj + l 



Cj-l 



Cn+j — 1 
Cn+j-2 



(2.5) 



If the conditions (12. 3p and (12. 4p are fulfilled, the polynomials Pn{z) satisfy the recurrence 
relation (see, e.g., [12]) 



Pn+l{z) + {dn - z)Pn{z) = zbnPn-l{z), n > 1, 

where the recurrence coefficients are 



A, 

Pn{0) 



T„,A 



/O, n = 0,l. 



n'-^n+l 



bn — dn 



hr. 



?n+lA„_i 



with T„ = An \ Note the important relation 

n = 1,2, 



n = 1,2, 



(2.6) 

(2.7) 
(2.8) 



hr 



A 



dn ^n— 1 



A2 



from which one can obtain expression for the normalization constant /i„ in terms of the recurrence 
parameters: 



hn = ll 



1=1 



di 



(2.9) 



There is a one-to-one correspondence between the moments c„, and the recurrence coeffi- 
cients bn-, dn (provided restrictions bndn / are fulfilled). 

We say that the LBP are regular if bndn 7^ 0. This condition is equivalent to the condition 

A„A(,i)y^O, n = 0,l,.... 

In the regular case there is a simple formula relating the biorthogonal partners Qn{z) with 
polynomials Pn{z) [12] : 



Qn{z) 



zPn+l{l/z)-Z^~^Pn{l/z) 
Pn{0) 



(2.10) 



In what follows we will use so-called rescaled LBP 
Pn{z)=q''Pniz/q), n = 0,l,... 
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with some non-zero parameter q. It is easily verified that the rescaled polynomials Pniz) are 
monic LBP satisfying the recurrence relation 

Pn+liz) + (dn - z)Pn{z) = zbnPn-l{z) 

with 

dn = qbn, bn = qbn- 

The rescaled LBP Pniz) differ from initial LBP Pn{z) only by a trivial rescaling of recurrence 
parameters. The moments c„ of the rescaled LBP are connected with initial moments c„ by the 
relation c„ = q^Cn- Note that the rescaled biorthogonal partners Qn{z) are transformed as 

Qn{z) = q-'^Qnizq). (2.11) 

There is a connection between the LBP and the restricted relativistic Toda chain [16]. Assume 
that LBP Pn{z;t) depend on an additional (so-called "time") parameter t. This mean that the 
recurrence coefficients bn{t), dn{t) become functions of the parameter t. We assume that the 
relation 

Pniz) = -^Pn-liz) 
dn 

holds for all n = 0, 1, ... . This ansatz leads to the following equations for the recurrence 
coefficients [16] 

dn = ¥^-:P^^ bn = bn(^- ^T^) ■ (2.12) 
"n+1 «n-l V"n "n-1/ 

For the corresponding moments c„(t) we have the relation 

c„ = c„_i, n = 0,±l,±2, 

Another possible ansatz [16] 

Pniz) = -bniPniz) - zPn-liz)) 

leads to the equations 

dn = -dnibn+i -bn), 6„ = -6„(5„+i - + - d„). (2.13) 

In this case we have for the moments the relation 

Cn = c„+i, n = 0,±l,±2, 

In spite of the apparent difference between equations (j2.12p and (j2.13p . it can be shown 
(see, e.g., [16]) that these two systems are both equivalent to the restricted relativistic Toda 
chain equations. The term "restricted" in this context means that it is assumed an additional 
condition 

bo = 0. 

This means that in formulas (|2.12p or (|2.13p we should assume n = 0,1, 2,.... For the non- 
restricted relativistic Toda chain equations ()2.12p or p.lSp are valid for all integer values of 
n = 0,±l,±2,.... 
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3 Laurent biorthogonal polynomials and QD-algorithm 

The (restricted) "discrete-time" relativistic Toda chain corresponds to the following ansatz for 
the moments 

Cn(t + h) = Cn+lit), n = 0,±1,±2, . . . , 

where h is an arbitrary parameter. We have the transformation formula for the corresponding 
Laurent biorthogonal polynomials 

Pn{z;t + h)=Pn{z;t)+bn{t)Pn^iiz;t) (3.1) 

and 

{dn - bn)Pn{z]t - h) = zPn{z;t) - Pn+l{z;t). (3.2) 

Formulas (j3.ip and (j3.2p can be interpreted as Christoffel and Geronimus transformations for 
LBP 

The corresponding recurrence coefficients are transformed as [32j 

dn{t + h) = dn~l ^J"^^ ^"^ , bn{t + h) = bn ^J"^^ (3.3) 
On — dn-1 On — fln-l 

(in r.h.s. of p.3p it is assumed the argument t for the coefficients dn)- These relations can 
be presented in a slightly different equivalent form as 

bndn = dn-lbn, bn — dn = bn+1 — dn, (3.4) 

where we have denoted 6„ = bn{t + h) etc for brevity. Relations (j3.4p describe so-called QD- 
algorithm for the two-point Fade approximation (see, e.g., [5] for details). In other words, 
the (restricted) discrete-time relativistic Toda chain is equivalent to the Q-D-algorithm for the 
two-point Fade approximation. 

Usually, this algorithm works as follows. We start from the given moments c„(t), n = 
0, ±1, ±2 . . . where the dependence on "time" is trivial: 

Cn{t + h) = Cn+l{t) 

and define the coefficient do(t) for all t = Iq + jh, j = 0, ±1, ±2 as 

co{t) 

The initial value to is not essential, usually it is assumed that to = 0, in this case we can write 
doit + jh)^dil^ = ^. 

Assume that bo{t) = for all t. Then at the first step we find bi{t) = b^^ for all t = jh from 



the second relation (j3.4p : 

"1 ~ "o "o 
Then we find dl from the first relation (|3.4p 



6 



S. Tsujimoto and A. Zhedanov 



(7) (7) (i) (i) 

This process can be continued to find 63 j ' • • • • "^^^ process is non-degenerate if bn dn / 



for all n and j. Then we obtain all sequences dii , h 



n 



0,1,2,... for j = 0,±1,±2,... . 



There is a remarkable connection with the Q-D-algorithm for the ordinary orthogonal poly- 
nomials [5]. Indeed, let us introduce the monic polynomials 



(3.5) 



where the polynomials Pn\z) are defined as PA^'iz) = Pn{z;hj). 
Then relations (13.11) and (13.21) become 



Wi^-'\z) =W(^\z) - fi^^wi^},{z) 



(3.6) 



and 



zWi^+^\z) 



(3.7) 



where 



fU) — kU+n-l) Jj) _ dj+n) _ Aj+n) 

J n "n 1 "n ""n 



These relations can be interpreted as Geronimus and Christoffel transforms for the orthogonal 
polynomials Wn\z). The compatibility condition between p. 60 and (13. 7p leads to the recurrence 
relation 

wi^Uz) +gU)wi^)(z)+u'^^wi^l,{z) = zWi^\z), 

which describes the three-term recurrence relation for the ordinary orthogonal polynomials 
Wn\z) where the recurrence coefficients are [5] 



Jn+l-: 



U) 



p{j) fiJ) 

J n • 



(i) (i) 

Moreover we have compatibility conditions for the coefficients , /A 



(i+i) f{j+i) ^ ^{j) f{j) 

J n 



^n-l Jn 



I f(j+l) 
-n ' J n 



(3. 



Relations (13. 8p coincide with those introduced by Rutishauser and describing the ordinary QD- 
algorithm [7]. It is easy to verify that relations (j3.8p are equivalent to relations (j3.4p for the 
two-point Q-D-algorithm. 

Thus starting from known solution Pn\z), bn\ dn^ of the discrete-time relativistic Toda 
chain (or, equivalently, two-point Q-D-algorithm) we can obtain a set of the ordinary orthogonal 

polynomials Wn\z) depending on additional "time" parameter j. Note that sometimes the 
introduced orthogonal polynomials Wn^iz) depending on an additional discrete parameter j are 
called the Hadamard polynomials [21 113)^1 . 

From the definition (13. 5p it follows that the orthogonal polynomials Wn\z) can be presented 
in determinantal form as 



rrU) 



Cj + l 
Cj+2 



Cj+2 
Cj+3 



1 z 



Cn+j+1 
Cj+n+2 

C2n+j 



^The authors are indebted to A. Magnus for drawing their attention to these references. 
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(7) 

where Hn stands for the Hankel determinant 



Cj+l Cj+2 
Cj+2 



Cj+n+1 



Cj-\-n . . . Cj-\-2n—l 

Clearly we have the relation 

H^^) = (-1)"("-i)/2a(,"+-''). 

Thus the orthogonal polynomials Wn\z) are orthogonal 

where the linear functional r(j') is defined by the moments 

r^^') = (r(^'),z") =c„+,+i, n = 0,1,2,..., j = 0, ±1, ±2, . . . . 

(7) 

The normalization constant qii has the expression 



H. 



(i) 



(j+n+l) 



TT\J 



It would be instructive to interpret p.ip and (j3.2p in terms of so-called bilinear technique 
by using the determinantal identities. This technique is standard in the theory of integrable 
systems. 

As a first step, we give a compressed expression to (i„, — 6„ as 



^(0) ^(0) ' 



which can be derived from the determinantal identity, or Jacobi identity, for the Toeplitz deter- 
minant: 

AHiiAEli = (A(f))^-A(r^)Ari). 
Then the relations (|3.1|) and (|3.2|) can be transformed to the following bilinear equations, 

(7) 

respectively, where the functions On are defined by 



a. 



U) 



S3) JJ) 



-0 



'1 



Ji) Ji) 
Ji) ^(i) 

1 Z 



.(i) 



(Note that is proportional to the Laurent biorthogonal polynomial Pn'{z). 



{3), 
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4 Laurent and Baxter biorthogonal polynomials 

There is an alternative (but essentially equivalent) approach to biorthogonal polynomials pro- 
posed by G. Baxter [1]. The pair Pn{z)^ Qn{z) of the biorthogonal polynomials is defined in 
this approach by means of initial conditions = Qq = 1 and the following recurrence system 

P„+i(z) = zP„(z)-eWg:(z), Q„+i(z) =zg„(z)-e(f)p„*(z), (4.1) 

where en '^^ are some complex coefficients. It is clear that = — P„_|_i(0), tn^ = — (5„+i(0). 
Notation P*{z) is standard for so-called reciprocal polynomials, i.e. P^{z) = z'^Pn{l/z), Q* (z) = 
z"'Qn{^/ z)- Assume that en^en\l — Sn^Cn^) 7^ (this is the nondegenerate case). Then, 
excluding Qniz) from the system (j4.ip we arrive at a 3-term recurrence relation for the polyno- 
mials Pn{z): 

Pn+liz) + dnPniz) = z{Pn{z) + 6„P„_i(z)) 

coinciding with (j2.6p . where 



(1) (1) 

Clearly, polynomials Qn{z) satisfy similar relations with interchanging superscripts 1, 2. 

Conversely, assume that we have the nondegenerate Laurent biorthogonal polynomials Pn{z) 
satisfying (12. 6p . We can construct their biorthogonal partners Qn{z) by (12.10p . Then it is 
elementary to verify that polynomials Pn{z), Qn{z) satisfy system (14. ip with = — P„+i(0), 
Cn = — Qn+i(0). Sometimes system (j4.ip is more convenient for analysis due to apparent 
symmetry between polynomials Pn{z), Qn{z) and corresponding coefficients en\ . Note 
also that the Laurent and Baxter biorthogonal polynomials in turn are equivalent to the so- 
called Laurent orthogonal polynomials proposed by Jones and Thron [15]. The Jones and Thron 
polynomials contains terms z'' with both positive and negative degree k. For details of this 
equivalence see, e.g., [12] and [22] . 

There is an important special case when all the Toeplitz determinants are positive A„ > 
and moreover the moments satisfy the condition 



(as usual, c„ means complex conjugation of c„). In this case the biorthogonal partners Qn{z) 
coincide with complex conjugated polynomials Qn{z) = Pn{z) and there exists nondecreasing 
function (t{9) of bounded variation on the unit circle such that the orthogonality relation 

Pn{e'')Pm{e-'')da{e) = hJnm (4.2) 



holds. I.e. in this case we have polynomials Pn{z) which are orthogonal on the unit circle 
(abbreviated as OPUC [23]). Historically, these polynomials were introduced first by Szego [29] 
and are called the Szego polynomials orthogonal on the unit circle. They satisfy the recurrence 
relation 

Pn+i{z) = zPniz) - a„z"P„(l/z), (4.3) 

where the coefficients = —Pn+i{0) are called the reflection (or Schur, or Verblunsky, ...) 
parameters. The relation (j4.3p was first derived by Szego himself [29]. The reflection parameters 
are complex numbers satisfying the important inequality 

|a„|<l, n = 0,1,2,.... (4.4) 
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In fact, condition (j4.4p is equivalent to the condition of positive definite Toeplitz forms A„ > 
or to existence of a positive measure on the unit circle providing orthogonality property (j4.2|) . 

If, additionally, all the moments are real, then they satisfy condition c-n = Cn- In this case the 
reflection parameters are real parameters satisfying the restriction —1 < a„ < 1, n = 0, 1, 2, . . . . 
The biorthogonal partners then coincide with initial polynomials Qn{z) = Pn{z)- It is easy 
to show that the measure da is symmetric with respect to real axis in this case, namely the 
function (t{9) satisfies the condition a{2TT — 9) + cr{9) = const. 

For further details concerning theory of OPUC see, e.g., |1H I24j. 

5 Frobenius elliptic determinant formula 
and biorthogonal functions 

Assume that Vi, Ui, i = 0,1, .. . are two arbitrary sequences of complex numbers. Let 

Hn = det ||(7ij||i,j=0,...,r!,-lj 

where 

a{ui + Vj + f3) 

9ij = —( \ wm exp(7i^^i + l2Vj), 

a{ui + Vj)a{P) 

where (t(z) is the standard Weierstrass sigma function (see, e.g., [H [30]) and (3, 71, 72 are 
arbitrary. 

Recall that the Weierstrass sigma function is defined by the infinite product [1] 



where the product is taken over all points of the lattice s = 2muji + 2m'uJ3, m, m! = 0, ±1, ±2, . . . 
excluding the point with m = m' = 0. 2lo\ and 2a;3 are the so-called primitive elliptic periods. 
It is convenient to introduce the third period 2^2 = — 2wi — 2a;3 [T]. The Weierstrass sigma 
function possess quasi-periodic properties [1] 

a{u + 2ijja) = — exp{2r]a{u + uJa))iy{u), a = 1,2, 3, 

where the constants ija are defined as 

110 = 0.^0), a = 1,2, 3 

and C{u) = a'{u)/a{u) is the Weierstrass zeta function [1]. 
We have 



(t{U + V + (3)Y[ cr{ui - Uj)a{vi - v 



Hn = , ( ^ ^ exp(7ii7 + 72^) (5.1) 

i,j 

n—l n—l 

(we denote U = Ui, V = ^ Vi for simplicity, moreover it is assumed that the upper limit 

i=Q i=0 

for i,j in the products is n — 1). 

Formula (j5.ip was obtained by Frobenius in [9]. A simple elementary method to derive formu- 
la (jS.ip can be found in [3]. Frobenius and Stickelberger derived also in [H] several other explicit 
formulas for "elliptic determinants" in connection with the theory of rational interpolation. 
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Let (pkix), ipki^), k = 0,1, . . . (we assume that 4>o = 'ipo = 1) be two sets of functions in some 
argument x. Assume that there exists a Unear functional C such that 

{C,cl)j{x)Tpi{x)) = Qij. 

The hnear functional C is defined on the space of functions constructed from bilinear combi- 
nations of the type 

f{x) = ^ Cik(t)i{x)ll}k{x) 
i,k=0 

with arbitrary coefficients Cik- 

Introduce the following functions 



Pnix) 



1 



and 



Qnix) 



1 



500 501 

510 511 

5n-l,0 5n-l,l 

0o(x) (t)i{x) 



500 510 

501 511 

50,n-l 5l,n-l 

ipoix) ^pi{x) 



50n 
5ln 

5n— l,n 
(pn{x) 



5n0 
5nl 

9n,n—l 
i'nix) 



(5.2) 



where 



An — Hn — det ||5ij||ij=0,...,n-l- 

By construction, these functions are biorthogonal 

1 Pn(,x')Q mi^xy) — ^n^nm 

with respect to the functional C, where the normalization coefficients hn are 

hr,. = — . 



(5.3) 



Expanding the determinant in (|5.2|) over the last row we have explicit expression for the poly- 
nomial Pn{x): 



n—k^ 



Pnk<Pk[x) 



k=0 



where 



Pnk 



Hn{k) 



The auxiliary determinants Hn{k) are defined by canceling the kth column, i.e. 

Hn{k) = det ||5rij(fc)||ij=o,...,n-l, 
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where 



9ijif^) 



a{ui+Vj{k) + p) 



a{ui + Vj{k))a{f3) 
Here the sequence Vi[k) is defined as 



exp(7iiii +-i2Vj{k)). 



Vi{k) 



Vi \i i < k, 
Vi+i if i > k. 



Thus the determinant Hn{k) is obtained from the determinant Hn by replacing sequence Vi with 
the sequence Vi{k). (By definition Hn{n) = Hn and Vi{n) = Vi.) But formula (|5.ip is valid for 
any sequences Uj, Vi. Hence we can calculate all the determinant Hn{k) explicitly. Omitting 
obvious calculations we present the result 



Pnk 



— p72(fn-ffe 



^a{U + V + Vn-Vk + P) 
a{U + V + l3) 



-|-|- a{ui + Vk) 
k\ fj^^ a{ui + Vn)' 



where 



n-l 

n ^(^" 

i=0 



n 

Vi) n cr(^i - ^fc) 
i=k+l 



are "generalized binomial coefficients" . Similar expression can be obtained for the biorthogonal 
partners Qn{x) if one replaces the parameters Vi with Uj. 

In case when the sequence vj is linear with respect to j: vj = wj+^ we obtain the conventional 
"elliptic binomial coefficients" |10| : 



[n]l 



[k]\[n-k]\ 



-1) 



-n k 



where [x] = a{wx)/a{w) is so-called "elliptic number" and [x]k = [x][x + 1] ■ • • [x + k — 1] 
is elliptic Pochhammer symbol. Note that usually the elliptic number is defined in terms of 
the theta function [x] = 6i{wx)/9i{w) [1^, but for our purposes these definitions are in fact 
equivalent. 

We thus constructed an explicit system of biorthogonal functions Qn{x) starting from 

the elliptic Frobenius determinant. This system can be further specified by a concrete choice of 
the basic functions (f)n{x), ipn{x) and the linear functional a. Note that the idea to construct 
explicit families of biorthogonal functions directly from corresponding Gram determinants is due 
to Wilson [31]. For general biorthogonal rational functions the determinant representation can 
be found e.g. in [25] and [6]. 



6 Laurent biorthogonal polynomials 
from the Frobenius determinant 

In what follows we will assume that the period 2uJi is a real while the period 2uj3 is purely imagi- 
nary. This means that the fundamental parallelogram is a rectangle. Such choice is standard 
for many practical purposes because in this case the function (t{x) takes real values on the real 
axis X [1]. This is important for existence of a positive orthogonality measure on the unit circle. 
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Put 

71 = 72 = 7> 



-iw + a, Vj = jw, 



where w is an arbitrary real parameter which is incommensurable with the real period 2uJi over 
the integers, i.e. we will assume that 

wNi / UJ1N2 (6.1) 

for any integers A^i, Then for the entries of the Frobenius matrix we have 

_ (7{w{j -i)+f3 + a) 

a{w{j-i)+a)a{P)- 

This matrix has the Toeplitz form. We can therefore define corresponding monic Laurent 
biorthogonal polynomials by the formula 



Pniz) = -- 



Co 



Cl 

Co 



C-n+1 C-n+2 
1 Z 



Cn.-l 



Cl 



(6.2) 



where the moments are defined as 

^ _ „ _ „^wn+^c (^{wn + /? + g) 

Cn — i/0,n — C , , \ f o\ 

a[wn + a)(T[p) 

and the Toeplitz determinant A„ is defined by (15. 3p . 

As in the previous section, define the elliptic numbers [x] as 

[x] = a{'wx)/cr{w), 
and the elliptic Pochhammer symbol 

[x]n = [x][x + I] ■ ■ ■ [x + n - 1]. 
The elliptic hyper geometric function is defined by the formula 



(6.3) 



r+lEr I r*; Z 



b' 



s=0 



[Qljs[Q2js • • • K+lJs Ms(s-l) 

mhiUh2]s---[br]s 



(6.4) 



where 



r+l 



i=l 



i=l 



We have 



Proposition 1. The Laurent biorthogonal polynomials defined by formulas ()6.2p and ()6.3p are 

expressed in terms of the elliptic hypergeometric function: 



Pn{z)=BnzE2 



; ze 



-7U1 



-n, a + 1, —[a + l)n — /? + 1 
a + 1 — n, — (q + l)n — /3 ' 

where a = aw~^ , [3 = (3w~^ and 

^ ^^^wn [-Q]n [an + /3 + n] 
[a + 1]„ [an + /3] 

is the coefficient to provide monicity Pn{z) = z^ + 0{z^~^) of the polynomials Pn{z). 



(6.5) 



(6.6) 
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Remark. The parameters of the elliptic hypergeometric function in our case satisfy condition 
1 + 61 + 62 = fli + 02 + as 



and hence M = in the definition of the hypergeometric function (j6.4p . Our definition of 
the elliptic hypergeometric function is in accordance with the conventional one [lUl \Tl\ . The 
main difference is replacing the theta functions with the Weierstrass sigma functions. This 
replacement leads to appearance of the additional factor e^^^i^-'^). Indeed, there is relation 
between these functions [1] 

a{z) = const • exp ( ) 6i{z / {2loi)) 



(the constant factor is not essential because it is canceled in all expressions for elliptic hy- 
pergeometric series). Using this relation we can replace all sigma functions with the theta 
functions 0i{z) which leads to formula (j6.4|) . 

Now we calculate the normalization coefficients /i„ directly from Frobenius formula ()5.ip : 

An fT(a) a{an + (3) [-a + l]n[a + 1]„ ' 

In what follows we will assume the following restriction a 7^ wm for any integers m. Indeed, 
otherwise the normalization coefficient hn becomes singular and we have a degeneration. 

We observe also that the determinant defined by (12. Sp with j = 1 is obtained from A^^^ 



by the shift of the parameter a ^ a + w because c„+i(a) = Cn{a + w). Thus in general we have 
the important formula 

A(f)(a) = Ania + jw). 

In particular, we have 

.(1) ^nli _ Tn+i _ e^("+") a{{a + w){n + l)+ (3) [n]!^ 

Tn a{a + w) a{{a + w)n + 0) [-a]„[a + 2]„' ^ ' ' 



Formulas ()6.7p and ()6.8p allow us to find explicit expressions for the recurrence coeffi- 
cients bn, dn- 

Indeed, from (12.71) and (12.81) we have 



d - - - n][/3 + (g + l){n + l)][/3 + an] 

hn ^ [a + (n + l)][/3 + (a + l)n][/3 + d(n + l)] 

and 

6„ = = ,-7^ [n]^[/3 + (d + l)(n + l)][/3 + a(n-l)] ^ ^^^^^^ 

K~i + {a + l)n][j3 + an\[a + n][a + n + I] 

We thus obtained a new explicit example of the Laurent biorthogonal polynomials which have 
both explicit expression in terms of the elliptic hypergeometric function (j6.5p and explicit re- 
currence coefficients (|6.9p . (jG.lOp . 

As a by-product, we have also obtained a new explicit solution of the discrete-time relativistic 
Toda chain or, equivalently, a new explicit solution of the two-point Q-D-algorithm. Indeed, the 
recurrence coefficients 6„, dn given by (16. 9p and (j6.10p provide an explicit elliptic solution of the 
two-point Q-D-algorithm (j3.4p with t = a, h = w. In turn, using correspondence (j3.5p we can 
obtain elliptic solution of the ordinary QL)-algorithm (13. Sp . or equivalently, the discrete-time 
Toda chain solutions. As far as we know these solutions are new. 

In order to find explicit (bi) orthogonality relation for these polynomials we need first the 
explicit Fourier expansion of the elliptic functions of the second kind. We will do this in the 
next section. 
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7 Fourier series of the elliptic functions of the second kind 

Assume that f{z) is the simplest ehiptic function of the second kind [T] 

^(£±^e- (7.1) 

a[z + a) 

with some complex parameters k, P, a, 7. The function f{z) is quasi-periodic with respect to 
periods 2a;i, 2ujs: 

f{z + 2uji)=fiif{z), f{z + 2uJ3) = fi3f{z), (7.2) 

where fii = g'^viP+^'^ii ^ = g2»?3/3+2t^37^ demand that function /(z) be purely periodic with 
respect to the (real) period 2wij: 

f{z + 2u;d) = f{z), 

where j = 1, 2, . . . is an arbitrary positive integer. This leads to the condition = 1 or 

j(a;i7 + r]i(3) = inm, (7.3) 

where m = 0, ±1, ±2, .... Note that for j > 1 we should avoid the values m = 0, ±j, ±2j, . . . 
because they correspond to pure 2LiJi-periodicity. Of course, it is assumed that m and j are 
coprime, i.e. //i is a primitive root of the unity of order j: 

^1 = e J . 

Moreover, we assume that a = — oq ~ ^ckij where both parameters ao,i are real and are 
restricted by conditions 

< ao < 2(Ji, 0<ai<2ja;3j. (7.4) 

Conditions (j7.4p mean that the parameter —a lies within the fundamental parallelogram (i.e. 
rectangle in our case). If a takes values beyond this parallelogram, it is possible to reduce it to 
canonical choice ()7.4p using shifts by periods 2a;i, 2a;3. Due to quasiperiodicity property of the 
function f{z) this will lead only to redefining of the parameter 7. Moreover we assume that the 
imaginary part — qi of a is nonzero. This assumption is very natural if we would like to avoid 
singularities of the function f{z) on whole real axis. Equivalently, one can present a in the form 

a = — ao — 2i/ll'3, (7.5) 

where < < 1 is a fixed parameter which describes the relative value of the imaginary part 
ai = —2ivoj'i with respect to the imaginary period 2a;3. 

Thus we have the function f{z) which is periodic and bounded on the whole real axis. It is 
possible therefore to present f{z) in terms of the Fourier series 

00 / ■ \ 

f{z)= V ^„exp ^ . 7.6 

Our problem now is to calculate the Fourier coefficients A^. 
By definition, 

1 „, s ( —2T\inz\ , „ , , 

A„ = - y /(z) exp ( j T = 2jui (7.7) 
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(the integral is well defined because by our assumptions the function /(z) has no singularities 
on the real axis). 

In order to calculate the integral in (j7.7|) we exploit standard method of contour integration 
(see, e.g., [1] for calculation of the Fourier expansion for Jacobi elliptic functions). Choose the 
contour F as the rectangle with vertices (0, 2ju;i, 2ju;i + 2a;3, 2103) (i.e. the horizontal length 
is 2ju)i and vertical length 2|a;3|). 

We have (the contour is traversed counterclockwise) 

where f-^, Jg correspond to horizontal sides of the rectangle, and integrals f^, correspond to 
vertical sides. 

Due to periodicity property f{z + 2juJi) = f{z) we have /2 + /4 = 0. For the two remaining 
horizontal integrals we have 

- 

Jl 

and 



f /-St^a+T /-2Trinz\ , 



Making the shift z ^ z + 2uj3 and using quasi-periodic property ()7.2p we have 



and thus 

/W/rexp (^Z^) = (1 - «exp (-i^)) A^. (7.8) 

Hence, in order to calculate the Fourier coefficient An we need to calculate the contour integral 
in l.h.s. of (j7.8p . This can be done by standard methods of residue theory. 

Indeed, inside the contour F the function /(2;)exp ( ~^™^ ) has only j simple poles located 
at points 

Zs = ao + iai + 2sa;i, s = 0, 1, . . . , j - 1. 
At zq = —a = ao + iai the function f{z) has the residue 

r = Ke-^"cj(/?). 
At Zs the function f{z) has the residue 

Ts = fJ^lr. 

Hence we have that the residue Rn of the function /(z)/Texp (^) inside the rectangle F will 
be 

i?„ = ^^^A^?e'^^^ = ^^(l + g + g2 + ... + gJ 1), 

s=0 
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where 



■Kin rp ( 2iri(m — n) 

X = ~- — ) Q = = exp 



If n 7^ TTimodj then Rn = 0. Nonzero value of the residue will be only for n = m + jt, 
t = 0, ±1, ±2, .... In this case 

= = ^exp f ^1 exp f 



Comparing with ()7.8p we get 
2'jriRn 

An = ^— ^, n = m,m±j,m±2j,... 

l-^3exp(-^) 

and 

An = 0, if n ^ m mod j. 
We can simplify this expression using the Legendre identity [1] 

ITT 

m^s - m^i = y 

which is valid if Im(u;3/cJi) > 0. Also we use the notation [1] 

n = exp . 

V '^i / 

In our case when loi > 0, ioo^, < we have that < /i < 1 (this is so-called normal case for the 
elliptic function [1]). 
We then have 



and 



Rn = Roexp[ '-)h ^ 

J^i J 



where 



Kfj(/3) / a[3rji 
Rq = exp ' 



2uJi \ uJi 

and we took into account relation (j7.5p . 

Thus for n = m + jk, k = 0, ±1, ±2, ... we have 

An = ^ (7-9) 

1 — e 

and An = if n 7^ m modj. 

Recall that j is a fixed positive integer - the order of the root of unity ^1, while m is a fixed 
nonnegative integer (lesser than j) coprime with j. Thus for large j the nonzero coefficients An 
are more rare then for small j. 
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There are two important simplest cases: 

(i) if J = 1 and m = 0. This case corresponds to the period 2uji. Then the Fourier coefficients 
An are nonzero for all ?i = 0, ±1, ±2, . . . and we have 

2TriRoexp (-^^) h-^'^" 
1 - e ^1 

(ii) if j = 2 and m = 1. This case corresponds to the period \lo\. In this case all even Fourier 
coefficients are zero = and for the odd Fourier coefficients we have 



27rii?oexp (-^^) /i-^""^ 
1 - e~~/i-2" 



^2n+l — -—p 



Note that in all cases the Fourier series ()7.6p converges inside the strip —v\ < lm{z) < V2, 
where 

vi = 2\u;3\{l - u), V2 = 2\uj3\iy. 

This results follows from standard theorems concerning asymptotic behavior of the Fourier 
coefficients An and for n ^ oo jl]. The parameters vi, V2 are positive as follows from the 
inequality < < 1. These conditions are very natural because the boundary lines Im(z) = 
2;/|w3| and Im(2;) = 2(i/ — l)|w3| of the strip pass through the poles of the function (j){z). Note 
that for u = 1/2 (i.e. when the pole of the function (j){z) lies on the horizontal line Im(z) = {losD 
we have the strip symmetric with respect to the real line: |Im(2;)| < |a;3|. The latter case 
correspond, e.g., to the Jacobi elliptic functions sn(z;/c), cn(z;A:), dn{z;k) [1]. 

8 Explicit biorthogonality relation 

In this section we obtain explicit biorthogonality property of the obtained Laurent biorthogonal 
polynomials. 

To do this we need to find explicit realization of the moments Cn given by formula (16. 3p . We 
note that 

Cn = f{wn), 

where f(z) is the elliptic function of the second kind ()7.ip (in our case k = l/cr(/3) but the 
constant k does not play any role in formulas for the polynomials Pn{z) and their recurrence 
coefficients 6„, 

Assume first that the parameter 7 is chosen to provide the periodicity of the function f{z) 
with period 2uJij, j = 1,2, . . . . Then we have the Fourier expansion (|7.6p from which one obtains 

A.expf— — j = (8-1) 

s=—oo ^ / 

where 

/ i-Ksw\ 

= exp , s = 0,±l,±2,... 8.2 

V J^i J 

is an infinite set of points belonging to the unit circle \zs\ = 1. Due to condition (|6.ip we have 
that all these points are distinct Zs zt \i t ^ s and hence they are dense on the unit circle. 
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Prom (jS.ip it follows that the moments Cn are expressible in terms of the Lebesgue integral 



1 

2^ 



27r 







over the unit circle j^l = 1, where ^{9) is a (complex) function of bounded variation on the 
interval [0, 27r] consisting only from discrete jumps Ag localized in the points Og given by (|8.2p . 

Thus we found explicit realization of the moments Cn and hence we immediately obtain 
biorthogonality relation for our Laurent biorthogonal polynomials 



AsPn{Zs)Qm{'^/Zs) = hn5nm, (8.3) 



where Qn{z) are biorthogonal partners ()2.2p with respect to polynomials Pn{z). The Fourier 
coefficients Ag play the role of discrete weights in this biorthogonality relation. Hence we have 
obtained 

Proposition 2. In the periodic case f{z + 2ujij) = f{z) the elliptic polynomials (j6.5p Pn{z) are 
biorthogonal (I8.3p on the unit circle \z\ = 1 with respect to a dense point measure with weights Ag 
given by expression 



Note that the biorthogonal partners Quiz) in our case can be found explicitly in terms of 
the elliptic hypergeometric function. Indeed, from (j2.2p we see that the polynomials Qn{z) are 
Laurent biorthogonal polynomials corresponding to the "reflected" moments c„ = c_„. From 
explicit expression (|6.3p it follows that the moments c—n are obtained from the moments Cn 
by reflection of the parameters a —a, jS — > — /?, 7 — > —7, whereas the parameter w remains 
unchanged (under such procedure we obtain the moments — c_„ but any constant common factor 
in front of moments leads to the same polynomials Qn{z))- Hence we can obtain expression 
for the polynomials Qn{z) from the expression (j6.5p for polynomials Pn{z) by reflection of 
parameters a, (3, 7: 



Qniz)=Rn3E2 



-n,l-a,{a-l)n + P + ^^^^ 
1 — n — a, {a — l)n + (3 ' 



where the coefficient Bn is obtained from corresponding coefficient (j6.6p by the same 
reflection of the parameters q, /?, 7. 

Thus both polynomials Pn{z) and their biorthogonal partners Qn{z) have similar expressions 
in terms of elliptic hypergeometric function. 

So far, we assumed that the function f[z) is periodic with the period 2u)ij. This assumption 
means that the parameter 7 should satisfy condition (j7.3p . Parameters a and /3 are assumed 
to be arbitrary (with the only condition (17. 4p ). What happens if the function f{z) is not 
periodic, i.e. if the parameter 7 is arbitrary? It appears that this general case can be easily 
reduced to the already considered. Indeed, assume that we change the parameter 7, i.e. assume 
that the parameters a and (3 remain the same but 7 = 7 + X, where x is arbitrary complex 
parameter. Then it is easily seen from explicit expression (16. Sp that the new Laurent biorthogonal 
polynomials Pn{z) are obtained by simple rescaling of the argument: 

Pn{z)=q^Pn{z/q), 

where q = e^^. This corresponds to transformation of the moments c„ = eq^Cn as seen directly 
from (|6.3p (the common constant e = e""' is inessential and can be put equal to 1). 
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Assume that we choose the parameter x such that the new function 
a[z + a) 

will be periodic with the period 2u;ij. This means that the parameter x should be chosen from 
condition (see (17. 3p ) 

^1(7 + X) = ^7rm/i> (8-4) 

where m is co-prime with j. 

Then the new polynomials Pn{z) will be biorthogonal on the unit circle according to above 
obtained proposition: 

00 

AsPniZs)Qmi'i-/Zs) = hndnm, (8.5) 

s=— oo 

where the spectral points Zg on the unit circle are given by (j8.2p and the weights As by (|7.9p . 
Note that the normalization coefficients /i„ remain unchanged under the rescaling transform as 
seen from (j2.9p . i.e. hn = hn- 

Taking into account that Qn{z) = q~^Qniz) (see (12. lip ) we obtain from (18. 5p the biorthogonal 
relation 

oo 

AsPn{Zs/q)Qm{q/Zs) = hnSnm- (8.6) 

s=— oo 

Relation (18. 6p means that for generic values of 7 polynomials Pn{z) and Qn{z) are biorthogonal 
on the non-unit circle \z\ = l/\q\ with respect to the same dense point measure. 

It is interesting to note that for every integer j = 1,2,... (i.e. for every period T = 2oJij) 
we can construct corresponding circle providing biorthogonality relation (j8.6p . Thus there exist 
infinitely many orthogonality circles for different values of the integer parameter j. 

For the radius r of the circle of biorthogonality we have from ()8.4p (recall that we assume 
parameter w to be real) 

r]i (3w 

r = ll\q\ = \e~\ le'^^l . 



9 Positivity of the measure and polynomials orthogonal 
on the unit circle 

Return to the case when the function f{z) is periodic with the period 2a;ij and consider an 
important special case when all the Fourier coefficients of the function f{z) are nonnegative 
An > 0. In this case all spectral points Zg belong to the unit circle = 1 and the measure on 
the unit circle is a positive nondecreasing function. 
We have < h < 1. Thus for n — > —00 we have 

An = 2mRoe~^ h-'^'^'' . 

It is seen that for positivity of An one should have 27rii?o = ^^Q^ where kq is a positive parameter, 
and for the real part of a we have the condition 



Oo = 2Joa;i, Jo = 0,±l,±2, 



(9.1) 
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Now for for n ^ oo we have 

An = -Koe^h'^^^-"^''. 
In this case we should have necessarily 

Re(/3) = (2Ji + Ji = 0, ±1, ±2, . . . . (9.2) 

It is easily seen that conditions (j9.ip and (|9.2p are also sufficient and so we have the 

Proposition 3. The Fourier coefficients are positive {up to inessential common factor) if and 
only if the real parts of parameters a, (3 satisfy conditions (I9.ip and (19. 2p . In this case the 
expression for the Fourier coefficients can he presented in the form 

^n = Ko— — 25", n = m + jk, k = 0,±l,±2,..., < v < 1, (9.3) 

-L ~r ^1 n 

TT Im(/3) 

and An = ifn^m mod (j), where ki = e "i is a positive parameter {as usual by Im(/3) we 
denote the imaginary part of (3). 

In this case we have positive dense point measure on the unit circle. It is well known that 
when the measure da is positive on the unit circle then biorthogonal polynomials become the 
orthogonal polynomials on the unit circle |29t [TTt I24j . In this case the moments c„ satisfy the 
restriction 



and moreover all the Toeplitz determinants are positive 

A„>0, n = l,2,.... 

The property c-n = Cn can be verified directly from the definition (|6.3p if the parameters q, 13 
satisfy conditions: 

a = 2Jia;i - 2i/a;3, j3 = {2Ji + l)oji + ij3i (9.4) 

(here (3i is an arbitrary real parameter). 

In this special case the obtained polynomials satisfy the Szego recurrence relation (j4.3p . The 
reflection parameters o„ are calculated as = — P„+i(0) and using already found explicit 
formula (j6.5p for polynomials Pn{z) we have = —Bn+i-, where Bn is given by (|6.6p (with q, (3 
satisfying restrictions ()9.4p ). From general theory it follows that in this case the reflection 
parameters should satisfy the restriction \an\ < 1- This property is not obvious from explicit 
expression for a„ in terms of elliptic Pochhammer symbols. 

If, in addition to positivity of An, we demand that the discrete measure should be symmetric 
with respect to the real axis we then obtain the condition j4_„ = An for all n = 0, 1, 2, ... . It 
is easily verified from explicit expression (j9.3p that this is possible only for j = 1 and j = 2. In 
the first case, when j = 1 the period T = 2lo\ and necessarily = 1/2 and k\ = 1, so that 

An = , . (9.5) 

K^ + h-'' ^ ' 

But the Fourier coefficients with expression (j9.5p correspond to the Jacobi elliptic function 
dn(z; /c) [30]. In this case the moments are c„, = dn('u;n; A;) and indeed satisfy the property 
C-n = Cn] the reflection parameters are very simple: a„ = An{w{n + 1); k) for the even n and 
ttn = —cn{w{n + l);k) for the odd n. 
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In the second case, i.e. when j = 2 we have the period T = Auji and necessarily = 1/2 and 

Ki = /i^"*^, so that 

A - ^0 

These Fourier coefficients correspond to the Jacobi elliptic function cx].{z;k) [30]. Again the 
moments Cn = CTi{wn; k) satisfy the desired property c_n = c„ and we have the polynomials 
orthogonal on the unit circle with simple reflection parameters: a„ = cn(w{n + l);k) for the 
even n and a„ = —dn{w{n + 1); k) for the odd n. 

These two explicit cases of OPUC with dense point measure were first considered in [33]. 
Now we see that there exists much wider class of explicit elliptic OPUC with positive dense 
measure on the unit circle. This class of OPUC contains essentially 3 arbitrary continuous 
parameters: w, Im(a) = — ai, Im(/3). We thus have two additional parameters with respect to 
the only parameter w in [33j . Note however, that if one demands that OPUC were real (i.e. they 
have real reflection parameters an and moments Cn) then nothing more general than "cn-" and 
"dn-" polynomials considered in [33] appear. 



10 "Classical" property of LBP 

Assume that Pn{z) are arbitrary Laurent polynomials satisfying 3-term RR 

Pn+l{z) + dnPn{z) = z{Pn{z) + 6„P„_i(z)) 

with some recurrence coefficients dn- 

For any sequence n = 0,1,2, .. . of complex numbers such that = we define the 
linear operator T> which acts on the space of polynomials in the argument z by the rule 

Then it is clear that the operator D sends any polynomial of degree n to a polynomial of degree 
n — 1 and moreover I'{1} = 0. In this sense the operator T) can be called as a generalized 
derivative operator. If /i„ = n then T> = dz coincides with the ordinary derivative operator with 
respect to the variable z. 

We say that LBP Pn{z) are P-classical if 

VPniz) = flnPn^liz), (10.1) 

where Pn{z) is another set of LBP satisfying the recurrence relation 

Pn+l{z) + dnPniz) = z{Pn{z) + 6„P„_i(z)) 

with some coefficients dn, bn and initial conditions Pq = 1, Pi = z — do. 

In |32j we considered the case of the ordinary classical LBP (i.e. with respect to the opera- 
tor dz) and derived necessary and sufficient conditions for existence of such polynomials. It 
appears that there exists many different types of such classical LBP. The simplest ones are the 
LBP constructed by Hendriksen and van Rossum [12|. The latter have explicit expression in 
terms of the Gauss hypergeometric function. 

Now return to our elliptic LBP Pn^z) = Pn{z;a,P,^,w) (we indicate dependence on param- 
eters a, P, 7, w for convenience) and consider the operator V with /i„ defined as 
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Then from explicit representation (j6.5p it is elementary to verify that the operator T) transforms 
these Laurent biorthogonal polynomials to the same family but with the sole parameter (5 
changed: 

VPniz;a,f3,-f,w) = finPn-iiz;a, f3 - a,-f,w) 

which means that our elliptic polynomials Pn{z;a, (3,^,w) are indeed "classical" polynomials 
with respect to the operator V. 

In particular, the choice a = —ws, (3 = uJi corresponds to the OPUC dn-polynomials [33] . 
Under action of the operator T> we obtain polynomials with a = — W3, (3 = loi + uj^ = —uj^- 
These polynomials correspond to the OPUC cn-polynomials [33] . Vice versa, action of the opera- 
tor T) on the cn-polynomials return them to dn-polynomials. In this case the coefficient /i„ is 
proportional to the Jacobi sn-function: = const • sn(tt;n; /c) [33] . 

One can repeat action of the operator T>. This leads to a chain of corresponding transforma- 
tions of polynomials -Pn(^): 

V^Pniz; a, (3, 7, w) = ^n/^n-l • • • fJ.n-m+lPn-m{z] a, (3 - ma, 7, w). 

One can consider the "/i-exponential" function E^{x) which is a formal solution of the opera- 
tor equation 

VE^{x) = E^{x). (10.3) 
Clearly we have a solution of the operator equation ()10.3I) in terms of the formal series 



In case of the elliptic fin (110. 2p we have 



s=Q 



Function (jl0.4p is closely related with so-called "theta analogue" of the exponential function 
proposed by Spiridonov in |26j^ . 

Obviously for an arbitrary complex parameter 7 we have 

We can also introduce "even" and "odd" /^-exponential functions which are ^-analogs of the 
hyperbolic "cosh" and "sinh" functions 

00 ^2s °o ^2s+l 



s=0 ^^^2 ■ ■ ■ /^2s ^ fJ'lfJ'2 ■ ■ ■ I-I-2S+1 

We have obvious relations 

C^{x) = (E^ix) + E^{-x))/2, S^{x) = (E^ix) - E^{-x))/2. 
These functions both have the same property 
VC^i^x) = 7'C^(x), V^S^i^x) = -i^S^{x) 



^The authors are indebted to V. Spiridonov for drawing their attention to this result. 
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with an arbitrary parameter 7, and hence for arbitrary parameters (3q, Pi the function f{x) = 
PoC^i'jx) + /?iS'^(7x) is a formal solution of the operator equation 

Note the obvious "intertwining" property of these functions: 

1^Cf,{jx) = jSf,{jx), VS^i-^x) = 7C^(7x). 

We can use these properties in order to construct formal generating functions for the cn and 
dn-circle polynomials. 

Indeed, let P^\z) and Prl^' be cn and dn-circle polynomials corresponding to the choices 
a = — ws, (3 = —UJ2 and a = — ^3, j3 = oji respectively. As shown in [33] these polynomials 
satisfy intertwining properties 

PPf )(Z) = PPr^(^) = /^nPfU^) 

and hence 

where one can choose fin = sn^wn) / sn{w) . 

Construct the generating functions for the polynomials P^\z) and Pn^\z) as the formal 
series 

F(C)^,.t) = y ^ ^" , F^^\z;t) = y * ^" ^'K (10.5) 

^ • • • /in ^ /il/i2 • • • 

We have obviously 

V,F^^\z;t) =tF^^\z;t), V,F^^\z;t) = tF^^\z;t), 
where notation means that the operator V acts only on the variable z. As a consequence 

VlF(<^^{z;t)=t^F^^\z;t), VlF^^\z;t) = t^F'^^\z;t). 

This property means that both functions F^^\z;t) and F^^\z]t) can be expressed in terms of 
"even" and "odd" /z-exponential functions with respect to the variable z: 

F^^\z;t) = ^o{t)C^{zt) + Ci{t)S^{zt), (10.6) 
F(^)(z;t) = 7]o{t)C^{zt)+7]i{t)S^{zt) (10.7) 

with some functions S,i{t), rji{t), i = 1,2. For these functions we have relations ??o(i) = ^1(^)1 
r]i{t) = ^o{t) which follow easily from intertwining relations. Remaining functions ^0(^)5 
can be found as follows. Put = 0. Then from definition (jlO.Sp we have 



where bP , B^n^ are corresponding normalization coefficients (16. 6p . On the other hand from 
(fToel) and (fTOTl) we have 

F^''\Q-t) = Ut), F(^)(0;t) = ei(t) 
because C/,(0) = 1, 5^,(0) = 0. 
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We thus have exphcit expressions for the functions £,Q{t), ^i(t) in terms of formal series 

m = E ^ ^1 w = E • 

^ /^l/^2 ■■■fJ-n ^ • • • /in 

As shown in [33] the coefficients -BI*^^ and si^^ are expressible in terms of the Jacobi elliptic 
functions as follows: 



^{C) _ _ (C) _ [ -cn(u;n), n even 
n n-1 I dn(?im), n odd 



and 



^(D) ^ _^{D) ^ ( -dn{wn), n even, 
" ""-"^ cn(u;7i), n odd, 

where ai'^'* and ai^'* are corresponding reflection parameters. 

We can thus present expressions for ^o(^) and ^i(t) in a more explicit form: 

y. cn(2wg)t^^sn^^(u') y. dn{w{2s + l))t^'+hn'^'+\w) 
sn(w)sn(2u;) • • • sn(27i;s) sn{w)sn{2w) ■ ■ ■ sn{w{2s + 1)) ' 

^ _ y. dn(2^;s)t^^sn^^(w) _ y cn{w{2s + l))t^'+hn'^'+\w) 
sn(w)sn(2u;) • • • sn(27i;s) sn{w)sn{2w) ■ ■ ■ sn{w{2s + 1)) 

Hence we were able to find explicitly the generating functions for the OPUC Pn^^ (z) and Pn^^ (z) 
in terms of /^-exponential functions (I10.6P , ()10.7p . The problem of convergence of these functions 
needs a separate investigation. 

11 Rational limit of the elliptic functions 

and corresponding Laurent biorthogonal polynomials 

In this section we consider the rational limit of the elliptic functions, when both periods 2u;i, 2lo2 
tend to infinity. In this case for the Weierstrass functions we have simple formulas [1] 

(^(z) = z, Ciz) = 1/z, p{z) = l/z^. 
Hence for the moments Cn (16. 3p we have (one can put w = 1 and 7 = without loss of generality) 

n + (5 + a 



(3{n + a) 



11.1 



In this case the elliptic Frobenius determinant A„ becomes the well-known rational Cauchy 
determinant [18] and for corresponding Laurent biorthogonal polynomials (j6.2p we have (cti = 
a + 1) 

/ \ R f-n,ai,-ain- [3 + 1 \ 

Pniz)=Bn3F2[ ^ ]Z] (11.2) 

\ ai — n, —a\n — p J 

with 

{—a)n an + f3 + n 



Br, 



{a + l)n an + (5 
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where 3F2{z) is the ordinary hyper geometric function [17] and (a)„ = a(a + 1) • • • (a + n — 1) 
is the Pochhammer symbol (shifted factorial). These polynomials satisfy the 3-term recurrence 
relation (12.61) with 



^ _ ia-n){l3 + ai{n + l)){f3 + an) 
" (ai +n)(/5 + ain)(/3 + a(n + l)) ^ 

and 

" {(3 + ain){f3 + an){a + n + l){a + n)' ^ 
Note that as in the elliptic case we should require that 
a / 0,±1,±2, ... 

because otherwise the normalization coefficients hn become singular at some n > which means 
a degeneration. 

As in the elliptic case the recurrence coefficients ()11.3p and (jll.4p provide an explicit so- 
lution of the two-point QD-algorithm (j3.4p with t = a, h = 1. Similarly, one can construct 
corresponding solutions of the ordinary Q-D-algorithm (13. 8p . 

The biorthogonal partners have the expression 

f-n,l-a,{a-l)n + (3 + l 
Qn{z) =Bn3,F2\ i\ , « 

In order to find the orthogonality measure for these polynomials we first note that the 
moments (jll.ip can be rewritten in the form 

Cn = /3-' + ^=r'+ci°), (11.5) 
n + a 

where the moments 

c^^^ = l/{n + a) = lim 

/3— >oo 

correspond to a special case of the "classical" Laurent biorthogonal polynomials considered by 
Hendriksen and Van Rossum [T2]. The moments c^f^ correspond to the recurrence coefficients 

a — n r? 

dn = ■ — , K = --, ■ w ■ — TT- 11-6 

a + n + 1 (a + n)(a + n + l) 

The polynomials Pn{z) are expressed in terms of the Gauss hypergeometric function [12] 

Pniz) = , . 2Fi[ ;z . 11.7) 

(a + 1)„ \a + l-n J 

Formulas (jll.7p and (jll.Gp are easily obtained from pi.3p . (|11.4p and (jll.2p in the limit /? — > oo. 
The biorthogonality property for polynomials (jll.7p was found in [12j 



Pn{z)Qm{'^/z)z°' ^dz = /l„5„m, 

C 



where the contour C is the unit circle and integrating path starts from z = where ^ = 1 
and stops at z = 1~ , where z'^~^ = e^'^*". 
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The moments (jll.Sp differ from the moments c„, only by the constant term (3~^. This means 
that corresponding orthogonahty measure on the unit circle has an additional concentrated mass 
at z = 1. Thus the biorthogonality relation for polynomials (|11.2p looks as 

Pn{z)Qmii/z)w{z)dz = hn6nm, 

where the weight function is 

Note that inserting the concentrated mass at point z = \ can be performed by the Geronimus 
transform of the classical Hendriksen-van Rossum polynomials (see [32] for details). 
The "classical" property (jlO.ip holds for the polynomials (jll.2p with 

n 

fJ-n = — ■ — ■ 
n + a 

In more details 

VPn{z;a,[3) = finPn-i{z; a, (3 - a), 

where Pn{z; a, P) stands for LBP (jll.2p with explicit dependence on the parameters a, (3. Note 
the operator "D in this case does not coincide with the ordinary derivative operator dz- Hence 
the polynomials (jll.2p provide one of the simplest nontrivial examples of the "classical" LBP 
with respect to a nonclassical "derivative" operator T). 

Consider also a special case when both a and (3 are purely imaginary parameters: 

a = isi, 13 = is2, siS2 > 0. 

Then we can put (obviously the moments Cn are defined up to a unnecessary common constant 
factor) by the expression 

_ n + i{si + S2) 
n + isi 

from which the condition c„ = C-n follows which means that the corresponding polynomials Pn{z) 
satisfy the Szego recurrence relation (j4.3p and one can expect that these polynomials are ortho- 
gonal on the unit circle with respect to a positive measure. According to the general theory this 
occurs if and only if the reflection parameters = — Pri+i(0) satisfy the restriction |an| < 1 for 
all n. 

from (jll.2p we see that = —Bn- Then 

|an-lP = BnBn- 

After simple calculations we get 

Tl? + $1 {Sin + S2Y 1 + 

where 



rin = n/si, f„ = n/{nsi + S2). 
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It is easily seen that the condition |a„| < 1 is equivalent to the condition siS2 > 0. In this case 
we can present the orthogonality property in the form 

=0 



where 



Using correspondence ()3.5p consider the ordinary orthogonal polynomials Wn\z) = Pn^^\z), 
where by Pn\z) we mean polynomials ()11.2p obtained under the substitution a — > a + j. 
We have 

(■) f-n,a + n + j + 1,1- p-n{a+j + n+l) \ 

V a + J + 1,-fJ - n[a + J + n + 1) J 

with 

^ _ {-a-j-n)n {a + n + j + l)n + (3 
{a + j + n + l)n (a + n + j)n + /3 

These polynomials are orthogonal with respect to the moment sequence 

" ^ ^n+a+j+1 

In order to get corresponding orthogonality measure we first consider the limit /3 — > 00. In this 
limit we have the polynomials 



V a+j+1 J 



with 

{-a-j-n)r 

r> — 



{a+j + n + l)n' 

These polynomials coincide with a special class of Jacobi polynomials which are orthogonal on 
the interval [0, 1] with the weight function 'w{x) = x"'^^: 







Indeed, for the moments we have 

1 1 

x°+%"da; = = lim ri^^. 

a + j + n + 1 fi^oo 

(i) — 1 

Hence, the polynomials Wn with nonzero /? correspond to adding a concentrated mass M = /? 

at the endpoint x = 1 of orthogonality interval for the Jacobi polynomials. I.e. the weight 

function for the polynomials Wn\x) with finite values of /3 has the expression 



wix) = x"-^^ + p-^5{x -1). 
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Such polynomials are called the Krall-Jacobi polynomials (see, e.g., |19]). These polyno- 
mials have a remarkable property: they are eigenfunctions of a fourth-order differential opera- 
tor [Ulig. 

Hence we can expect that our Laurent biorthogonal polynomials will satisfy corresponding 
fourth-order differential equation too. We first consider a more elementary case (3 = oo. Then 
the Jacobi polynomials WrP {x) satisfy the differential equation |17j 

x{l - x) + -(a + j + {a + j + 2)(1 - 2x)) 

+ n{n + a + j + l)W^^\x) = Q. (11.9) 

In order to return to the Laurent biorthogonal polynomials we need to shift the parameter 
J — > J — n. Then equation (IlLDh becomes 

z{l-z) ' +-{a + j-n + {a + j-n + 2){l-2z))- ^' 



dz'^ ' 2^ ' ' ^ ' ' " dz 

+ n{a + 3 + l)P^^\z) = {). 

We can rewrite this equation in the form of generalized eigenvalue problem 

AP^^\z) = KBPi^\z) (11.10) 

for two differential operators 

A = z{l- z)dl + (a + j + 1 - (a + j + 2)z)d,, B = {I - z)d, - a - j - 1, 

where A„ = n is the generalized eigenvalue. 

Consider now general case of finite values of the parameter p. The orthogonal polyno- 
mials (jll.Sp satisfy the differential equation of the 4th order [19\ [20] 

LWif\x) = XnW.jf\x), (11.11) 
where the differential operator L is 

L = (x(l - x)fd^ + x(l - x)(eix + eo)5^ + (1 - x){r]ix + r]o)dl + {dx + Co)5, 
with the coefficients 

^1 = -2(a + j + 4), ^o = 2(a + j + 2). 

Til = -{a + j + 2){a + j + 7)- 2/3, no = {a + j + 1)(q +j + 2), 

Ci = 2{a + j + 2){a + P + j + 1), Co = -2(a + j + l)(a + /3 + j + 2). 

The eigenvalue is 

An = ^{11 — l)(n — 2){n — 3) — ^in(n — l)(n — 2) — r]in(n — 1) -|- d^ 
= n{n + a + j + l){n'^ + n{a + j + 1) + 2/3). 

If we now substitute j j — n into above formulas, we return to the Laurent biorthogonal 
polynomials Pn^ {z) which satisfy a differential equation of the 4th degree which can be presented 
in the form 

(n^La + nLi + Lo)P^^\z) = 0, (11.12) 
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where Lq is a differential operator of the 4th order, Li of the 3-rd order and L2 of the second 
order. The operators Lq, Li, L2 do not depend on the parameter n. We thus see that the 
Laurent biorthogonal polynomials Pn\z) given by (jll.2p . satisfy the quadratic operator pen- 
cil equation (jll.l2p (with respect to the "eigenvalue" parameter n). When /? = 00 we have 
generalized eigenvalue problem (jll.lOp which is equivalent to a linear operator pencil. 

Hence the Laurent biorthogonal polynomials (|11.12p can be considered as biorthogonal ana- 
logs of the Krall-Jacobi orthogonal polynomials. They possess many useful properties of the 
Krall-Jacobi polynomials including the 4th order differential equation they satisfy. The main 
difference, however, is that in the biorthogonal case we have quadratic operator pencil equa- 
tion (I11.12P instead of usual eigenvalue problem (jll.lip . 

We thus see that many nontrivial properties of the elliptic Laurent biorthogonal polyno- 
mials are manifested already in the rational limit. Loosely speaking, one can say that elliptic 
biorthogonal polynomials (j6.5p are elliptic analogs of the Krall-Jacobi polynomials. 

Acknowledgements 

The authors thank L. Golinskii, A. Kirillov, C. Krattenthaler, A. Magnus, M. Rahman and 
V. Spiridonov for discussion. The authors are also indebted to the referees for careful reading 
the manuscript and many suggestions leading to improving of the text. 

References 

[1] Akhiezer N.I., Elements of the theory of elliptic functions, 2nd ed., Nauka, Moscow, 1970 (English transL: 
Translations of Mathematical Monographs, Vol. 79, American Mathematical Society, Providence, R.I., 1990). 

[2] AUouche H., Cuyt A., Reliable pole detection using a deflated gd-algorithm: when Bernoulli, Hadamard 
and Rutishauser cooperate, Numer. Math., to appear. 

[3] Amdeberhan T., A determinant of the Chudnovskys generalizing the elliptic Frobenius-Stickelberger- 
Cauchy determinantal identity. Electron. J. Combin. 7 (2000), Note 6, 3 pages. 

[4] Baxter G., Polynomials defined by a difference system, J. Math. Anal. Appl. 2 (1961), 223-263. 

[5] de Andrade X.L., McCabe J.H., On the two point Pade table for a distribution. Rocky Mountain J. Math 
33 (2003), 545-566. 

[6] Derevyagin M., Zhedanov A., An operator approach to multipoint Pade approximations, J. Approx. Theory 
157 (2009), 70-88. larXi^: 0802.3432 ', 

[7] Faddeev D.K., Faddeeva V.N., Computational methods of linear algebra, W.H. Freeman and Co., San 
Francisco - London, 1963. 

[8] Frobenius C, Stickelberger L., Uber die Addition und Multiplication der elliptischen Functionen, J. Reine 
Angew. Math. 88 (1880), 146-184 (reprinted in Ferdinand Ceorg Frobenius Gesammelte Abhandlungen, 
Vol. 1, Editor J.-P. Serre, Springer, Berlin, 1968, 612-650). 

[9] Frobenius C, Uber die elliptischen Functionen zweiter Art, J. Reine Angew. Math. 93 (1882), 53-68 
(reprinted in Ferdinand Georg Frobenius Gesammelte Abhandlungen, Vol. 2, Editor J.-P. Serre, Springer, 
Berlin, 1968, 81-96). 

[10] Gasper G., Rahman M., Basic hypergeometric series, 2nd ed.. Encyclopedia of Mathematics and its Appli- 
cations, Vol. 96, Cambridge University Press, Cambridge, 2004. 

[11] Geronimus Ya.L., Polynomials orthogonal on a circle and their applications, Amer. Math. Soc. Translation 
1954 (1954), no. 104. 

[12] Hendriksen E., van Rossum H., Orthogonal Laurent polynomials, Nederl. Akad. Wetensch. Indag. Math. 48 
(1986), 17-36. 

[13] Henrici P., Applied and computational complex analysis, John Wiley & Sons, New York, 1974. 

[14] Ismail M.E.H., Masson D.R., Generalized orthogonality and continued fractions, J. Approx. Theory 83 
(1995), 1-40, ,math.CA/9407213| 



30 



S. Tsujimoto and A. Zhedanov 



[15] Jones W.B., Thron W.J., Survey of continued fraction methods of solving moment problems, in Analytic 
Theory of Continued Fractions (Loen, 1981), Lecture Notes m Mathematics, Vol. 932, Springer, Berlin - 
New York, 1982, 4-37. 

[16] Kharchev S., Mironov A., Zhedanov A., Faces of relativistic Toda chain, Internat. J. Modem Phys. A 12 
(1997), 2675-2724, |hep-th/9606144| 

[17] Koekoek R., Swarttouw R.F., The Askey scheme of hypergeometric orthogonal polynomials and its q- 
analogue, Report 94-05, Faculty of Technical Mathematics and Informatics, Delft University of Technology, 
1994, |math.CA/9602214[ 



[18] Krattenthaler C, Advanced determinant calculus, Sem. Lothar, Combin. (1999), Art. B42q, 67 pages, 
,math.CO/9902004, 

[19] Littlejohn L.L., The Krall polynomials: a new class of orthogonal polynomials, Quaestiones Math. 5 (1982), 
255-265. 

[20] Littlejohn L.L., On the classification of differential equations having orthogonal polynomial solutions, Ann. 
Mat. Pura Appl. (4) 138 (1984), 35-53. 

[21] Magnus A. P., Special nonuniform lattice (snul) orthogonal polynomials on discrete dense sets of points, 
J. Comp. Appl. Math. 65 (1995), 253-265. 

[22] Pastro P.L, Orthogonal polynomials and some g-beta integrals of Ramanujan, J. Math. Anal. Appl. 112 
(1985), 517-540. 

[23] Ruijsenaars S.N.M., Relativistic Toda systems. Comm. Math. Phys. 133 (1990), 217-247. 

[24] Simon B., Orthogonal polynomials on the unit circle, American Mathematical Society Colloquium Publica- 
tions, Vol. 51, American Mathematical Society, Providence, R.I., 2005. 

[25] Spiridonov V.P., Zhedanov A. A., To the theory of biorthogonal rational functions, RIMS Kokyuroku (2003), 
no. 1302, 172-192. 

[26] Spiridonov V.P., An elliptic incarnation of the Bailey chain. Int. Math. Res. Not. 2002 (2002), no. 37, 
1945-1977. 

[27] Spiridonov V.P., Essays on the theory of elliptic hypergeometric functions, Russ. Math. Surv. 63 (2008), 
405-472, larXiv:0805.3"T35] 

[28] Suris Yu.B., A discrete-time relativistic Toda lattice, J. Phys. A: Math. Cen. 29 (1996), 451-465, 
so lv-int/9510007, 

[29] Szego G., Orthogonal polynomials, 4th ed., American Mathematical Society Colloquium Publications, Vol. 23, 
American Mathematical Society, Providence, R.I., 1975. 

[30] Whittaker E.T., Watson G.N., A course of modern analysis. An introduction to the general theory of infinite 
processes and of analytic functions; with an account of the principal transcendental functions, reprint of 4th 
ed. (1927), Cambridge Mathematical Library, Cambridge University Press, Cambridge, 1996. 

[31] Wilson J. A., Orthogonal functions from Gram determinants, SIAM J. Math. Anal. 22 (1991), 1147-1155. 

[32] Zhedanov A., The "classical" Laurent biorthogonal polynomials, J. Comput. Appl. Math. 98 (1998), 121- 
147. 

[33] Zhedanov A., Elliptic polynomials orthogonal on the unit circle with a dense point spectrum, Ramanujan J., 
to appear, larXiv:07l 1.4696 1 



